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Abstract 

Triple differential dijet cross sections in e ± p interactions are presented in the region 
of photon virtualities 2 < Q 2 < 80GeV 2 , inelasticities 0.1 < y < 0.85, jet transverse 
energies E^ x > 7GeV, E^, 2 > 5GeV, and pseudorapidities —2.5 < 77*, 772 < 0. The 
measurements are made in the 7*p centre-of-mass frame, using an integrated luminosity 
of 57 pb _1 . The data are compared with NLO QCD calculations and LO Monte Carlo 
programs with and without a resolved virtual photon contribution. NLO QCD calculations 
fail to describe the region of low Q 2 and low jet transverse energies, in contrast to a LO 
Monte Carlo generator which includes direct and resolved photon interactions with both 
transversely and longitudinally polarised photons. Initial and final state parton showers are 
tested as a mechanism for including higher order QCD effects in low Et jet production. 
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1 Introduction 



Jet cross sections in electron-proton collisions are successfully described by next-to-leading 
order (NLO) QCD calculations in most of the HERA kinematic range [1-6]. However, regions 
of phase space have previously been observed for which NLO predictions do not reproduce the 
data satisfactorily and leading order (LO) Monte Carlo simulations with different approaches to 
modelling higher order QCD effects are often more successful [3-8]. 

At HERA, a photon coupling to the incoming electron interacts with a parton from the 
proton. The measurement of dijet production is particularly suitable for the investigation of 
effects related to photon structure, which have been studied during the last two decades in 
e + e~ and ep collisions [9, 10]. In the "photoproduction" region of ep interactions, i.e. for 
Q 2 <C Aq CD , the interaction can be described by the sum of two contributions. In the direct 
photon process, the photon interacts as a whole with a parton from the proton, whereas in the 
resolved photon process, it behaves as a source of partons, which interact with partons from the 
proton. 

As the virtuality of the photon increases, the role of photon structure gradually changes. 
Whereas for quasi-real photons it is an indispensable theoretical tool, for Q 2 3> Aq CD the con- 
cept of a resolved photon is usually discarded and the data are analysed within the framework of 
perturbative calculations of direct photon processes. However, it has been argued [11-14] that 
the concept of the resolved photon is very useful phenomenologically for arbitrary Q 2 , provided 
the photon virtuality remains much smaller than some measure of the hardness of the process in 
which the photon participates. In our case this is satisfied for Q 2 <C E\, where E T denotes the 
jet transverse energy. Experimental evidence for the resolved virtual photon contribution has 
been found in a number of publications [7, 15, 16]. 

In this paper, new data on dijet production, obtained with the HI detector in the kinematic 
region of low to moderate photon virtualities 2 < Q 2 < 80 GeV 2 , are presented as triple 
differential distributions in Q 2 , the inelasticity y and variables characterising the final state jets. 
The data are compared with predictions within several theoretical approaches, differing in the 
way QCD effects are taken into account beyond LO, in order to identify which of them are 
successful in which regions. In doing that we use NLO calculations (i.e. including terms up 
to order act 2 ) as well as LO calculations supplemented with parton showers, which take into 
account leading logarithmic contributions to all orders. The effects of resolved virtual photons 
are studied for both transverse and longitudinal photon polarisations. 

The paper is organised as follows. After a review of various theoretical approaches to the 
description of interactions of virtual photons in Section 2, a brief description of the detector is 
given in Section 3. The data sample and event selection are specified in Section 4. Background 
subtractions, detector corrections and estimates of the measurement uncertainties are discussed 
in Section 5. The results are presented and discussed in Section 6. 
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(a) (b) (c) (d) (e) 



Figure 1 : Examples of diagrams for the production of at least two jets (incoming electrons and 
protons not shown): LO (a) and NLO (b) direct photon interactions; LO resolved photon inter- 
actions involving a quark (c) or a gluon (d) from the photon; NLO resolved photon process (e). 

2 Dynamics of Hard Processes in ep Collisions 
2.1 The DGLAP Approach 

The DGLAP approach uses parton distribution functions (PDFs) of the proton extracted from 
analyses of data on hard scattering processes. These PDFs depend on the factorisation scale fif, 
satisfy the DGLAP [17] evolution equations, and are sometimes called "integrated" to empha- 
sise the fact that their definition involves an integral over the virtualities of partons up to fij. 

In this approach, the cross section for dijet production in our kinematic region is given by 
the direct photon contribution, illustrated in Fig. 1 a,b, and expressed as 

a Dm ~ D Hp ® a *i > a °i = c *i a > + ° { S a2 s + ■ ■ ■ > (!) 

3 

where a e j denotes the cross section for a collision between the incoming electron and a parton j 
from the proton, , cf^, ■ ■ ■ are coefficients of an expansion of a e j in powers of a s and Dj/ p 
denotes the PDF of the proton. The term c^a s defines the LO cross section, whilst c^a s + 
cf}a 2 defines the NLO cross section. 

Recent analyses [3, 7, 15] of dijet cross sections in the region Aq CD < Q 2 < E\ have 
convincingly shown that the LO direct photon contribution lies significantly below the data. 
The NLO calculations, involving diagrams such as that shown in Fig. 1 b, bring the theoretical 
prediction closer to the data [3]. A recent HI analysis [4] indicates that even the NLO calcula- 
tions do not completely describe inclusive jet production at low Q 2 in part of the phase space. 
Large values of the NLO corrections, i.e. the ratio of NLO to LO predictions for the cross 
sections, and high sensitivity of the predicted jet cross sections to variations of the factorisation 
and renormalisation scales, strongly suggest the need for higher order (i.e. cfj, , • • • ) terms 
in Eq. (1). In the absence of a full calculation beyond NLO, some approximate procedure for 
resummation of the dominant higher order terms in Eq. (1) can be constructed. This procedure 
is based on the fact that in part of the phase space, the upper vertex of the diagram in Fig. 1 b can 
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be viewed as the splitting of the photon into a qq pair. Taking into account subsequent emissions 
of partons from this qq pair, these terms can be resummed into the PDF of the photon 1 , A/7* > 
as is done for instance in [11], and A/7£ m [18], where 7^ denotes the transversely and the 
longitudinally polarised virtual photon [19, 20]. 

Consequently, one can calculate the resolved photon contribution to the dijet cross section, 
corresponding to the graphs shown in Fig. 1 c, d, e, as 



RES 

(J 



fk®Y. Di hi 

k=T,L i,j 



Gij = cffa 2 + c\fa 3 s + 



(2) 



where i, j run over all partons in the photon and proton respectively, is the partonic cross 
section, cf^ at 2 defines the LO resolved photon cross section, c[fa 2 + c\fa^ the NLO resolved 
photon cross section and f T , fi denote the fluxes of transversely and longitudinally polarised 
virtual photons, respectively: 



f T (y,Q 2 ) 
h(y,Q 2 ) 



a 
27 
a 
27 



2(1 -y) + y 2 1 2m 2 e y 

y ~ Q 2 Q A 
2(1 -y) 1 
y Q 2 



(3) 
(4) 



The final dijet cross section is then given by the sum 2 of a D1R and <r RES . 



There is an interesting connection between the direct and resolved photon contributions. 
In a large part of the phase space, the NLO direct calculations (Fig. 1 b) can be reasonably 
well approximated by the sum of the LO direct (Fig. 1 a) and LO resolved (Fig. 1 c) photon 
contributions, provided the simplest expression, namely that given by the pure QED splitting 
of the photon into a qq pair, is used for the photon PDF [21]. In our kinematic region and for 
quark masses m 2 <^ Q 2 , the pure QED photon PDFs have the form 



(^75 Of-* Et) — 



a 

27 

0. 



3c^4x^(l 



(5) 

(6) 
(7) 



In Eqs. (5)-(7), denotes the electric charge of the quark q { and x 7 denotes the four-momentum 
fraction of the photon carried by the quark. The full expressions for the distribution functions 
from Eqs. (5)-(7), containing the exact Q 2 dependence with the correct threshold behaviour for 
Q 2 /m 2 — > 0, can be found in [14]. 

'This resummation actually yields the point-like (sometimes called "anomalous") parts of the photon PDF, 
which dominate in the kinematic region studied in this paper, where the hadronic (sometimes called "VMD") parts 
of the photon PDF are negligible. 

2 Care must be taken when adding the contribution of the LO resolved photon diagram (Fig. 1 c) to the NLO 
direct photon term (Fig. 1 b). To avoid double counting, the so called photon splitting term must be subtracted 
from the NLO direct photon contribution of the diagram in Fig. 1 b. In this paper the term "direct contribution" 
denotes the direct photon contribution before the subtraction of the splitting term. 
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2.2 The CCFM Approach 



The CCFM [22] approach uses the more general concept of an "unintegrated" PDF of the proton 
in the region of small Bjorken-a;. The virtualities and transverse momenta of the propagating 
partons are no longer ordered, as is the case for DGLAP evolution. Instead, an angular ordering 
of emissions is introduced in order to correctly treat gluon coherence effects [22]. Similarly to 
the case of the DGLAP scheme in Eq. (1), the cross section can be factorised into a partonic 
cross section and universal parton distribution functions according to [23] 

kT FACTORISATION f^ Z j2 u ~ ( X t.2 \ a /■ t.2 ~2\ /on 

J z z 

3 

where the partonic cross sections a e j have to be taken off-shell (i.e. dependent on the parton 
transverse momentum, k T ), (if is the factorisation scale related to the maximum angle allowed 
in the evolution, and the unintegrated parton distributions, Aj/ p {x, (tj), depend on an addi- 
tional variable, the k T of parton j. 

The CCFM evolution scheme provides a framework for the implementation of A; T -unordered 
initial state QCD cascades. The partons with the largest k T may come from any emission in the 
proton cascade, not necessarily from the hard subprocess as in the DGLAP framework. This 
can lead to events which have a similar topology to that of the resolved photon interaction in 
the DGLAP approximation [24], where hard partons are accompanied by softer partons from 
the photon remnant. 

The mean value of the proton momentum fraction x p , appearing as an argument of the 
unintegrated PDFs, is (x p ) ~ 0.03 in our kinematic region. Even though this value may not be 
small enough for the CCFM approach to be superior to that based on the standard integrated 
PDF and DGLAP evolution equations, it is interesting to compare the CCFM predictions with 
the data. Recently, such comparisons became possible using the CASCADE Monte Carlo (MC) 
generator. 



2.3 Programs for Dijet Calculations 

Several MC programs can be used for predictions of dijet cross sections, as discussed below. 

HERWIG [25] is a general purpose LO MC event generator, applicable to a wide range of 
hard processes and collisions, including direct and resolved photon interactions in the 
region of moderate Q 2 . It is based on LO cross sections in the DGLAP approach, inter- 
faced to leading log parton showers. Hadronisation is done via the decay of colourless 
clusters, formed during the hard scattering and parton shower stages. HERWIG is also 
able to model additional soft remnant-remnant interactions (the "soft underlying event"), 
accompanying the hard scattering process. The probability that a resolved photon event 
contains soft underlying activity has been adjusted so that the energy flow in and around 
the jets is well described. 

Only transversely polarised photons are included for resolved photon interactions in the 
default version of HERWIG. To investigate the contributions of resolved longitudinal pho- 
tons, we have modified HERWIG by adding the option of simulating the flux in Eq. (4). 
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Similarly the QED PDFs of the photon from Eqs. (5-7) have been implemented in HER- 
WIG in order to study the differences between the results obtained with the pure QED and 
the QCD-improved photon PDFs. The scales [i r and [if are set to a combination of Man- 
delstam variables [25], which roughly corresponds to l.lp T , where pr denotes transverse 
momentum of the outgoing parton from the hard interaction. 

RAPGAP [26] combines standard LO hard scattering matrix elements in the DGLAP ap- 
proach with parton showers and LUND string fragmentation, using JETSET [27,28]. 
Only the transverse virtual photon is considered in resolved photon processes. Soft un- 
derlying interactions are not modelled. 

DISENT [29] is a NLO DGLAP program for calculating dijet cross sections at the parton 
level. It is based on the dipole subtraction method [30] to regularise soft and collinear 
divergences. The factorisation scale, [if, was set to (E T ) = 9 GeV, since the program 
does not allow the user to set /x/ = Et for each point in phase space. In our kinematic 
region, the difference between the results obtained with [if — E T and [if = (E T ), tested 
using JETVIP, is very small. DISENT does not include resolved photon interactions. 

JETVIP [31, 32] is a NLO DGLAP parton level program which calculates both direct and 
resolved photon contributions. It is based on the phase space slicing method. We have 
performed systematic investigations of the stability of JETVIP calculations with respect 
to variations of the slicing parameter y c [21]. The direct contribution in our kinematic 
region is independent of y c to within 5% over the recommended range of its values 10~ 4 < 
He < 1CT 2 . The situation changes in the case of the NLO resolved photon contribution, for 
which the dependence on the y c parameter is significantly larger. The sum of NLO direct 
and NLO resolved JETVIP predictions varies by 30% in some bins for the recommended 
range of y c . We set y c = 0.003 in all JETVIP calculations, since the predictions are most 
stable around this value. 

We have observed non-negligible differences between the differential dijet cross section 
obtained with DISENT and the direct contribution from JETVIP (see Section 6.1). We 
have checked that this discrepancy is not caused by different input parameters or kine- 
matic cuts, as the leading order, i.e. 0(aa s ), contributions agree perfectly. We therefore 
test both programs in our analysis. 

CASCADE [33,34] is a LO MC event generator. It uses unintegrated gluon distribution func- 
tions of the proton, satisfying the CCFM equation, and correspondingly produces a k T 
unordered initial state parton shower. The LUND string model is used for fragmentation. 

The parameter settings of all MC programs and NLO calculations are summarised in Ta- 
bid. 

In the NLO calculations, JETVIP and DISENT, the massless partons entering the hard pro- 
cess are taken to be exactly collinear with the beam particles. On the other hand, the Monte 
Carlo generators HERWIG, RAPGAP and CASCADE generate initial state QCD parton show- 
ers, which influence the four-momenta of partons entering the hard process by generating the 
appropriate transverse momentum. Both initial and final state parton showers also provide more 
partons in the final state. These effects tend to produce more low Et jets in the LO models than 
in the NLO calculations. 



8 



Parameters HERWIG RAPGAP CASCADE DISENT JETVIP 



Version 


6.4 


2.8 


1.2 


— 


2.1 


Proton PDF 


CTEQ5L [35] 


CTEQ5L 


J2003 (set 1) [36] 


CTEQ6M 


CTEQ6M 


Photon PDF 


SAS1D [11]; 


SAS1D 


— 


— 


SAS1D 




[18] for 72 










Formula for a s 


one-loop 


one-loop 


one-loop 


two-loop 


two-loop 


Active flavours 


5 


5 


4 


5 


5 


PRSOF 


10% 










fir 


~ l.lpT 


Vipl + m 2 q ) 


V(Pt + ™ 2 q ) 


771* 


TP* 




~ l.lpT 


V(Pt + m V> 


given by 


9GeV 


TP* 








ang. ordering 






Hadronisation 


Cluster model 


LUND string 


LUND string 






mechanism 




fragmentation 


fragmentation 







Table 1: Parameters of the MC programs. The variable p T denotes the transverse momentum 
of the parton with mass m q outgoing from the hard interaction and E^ 1 is the energy of the 
jet with the highest transverse energy. The parameter PRSOF specifies the fraction of resolved 
photon events with soft underlying activity. 



2.4 Hadronisation Corrections 

The MC event generators have a clear advantage over the parton level calculations in incor- 
porating hadronisation effects. In order to estimate the hadronisation corrections to the NLO 
calculations, we use two different Monte Carlo models, HERWIG and LEPTO [37], and divide 
the cross sections obtained from these models for the complete hadronic final state by the cross 
sections predicted from the partonic final state after the initial and final state QCD parton show- 
ers. The hadronisation corrections determined by HERWIG also include corrections for the soft 
underlying event. The average values of the corrections obtained with HERWIG and LEPTO are 
applied to the NLO calculations as bin-by-bin correction factors and half the difference between 
the corrections obtained with the two models is taken as a hadronisation uncertainty in the NLO 
predictions. The hadronisation effects usually do not change the NLO predictions by more than 
5%, with the exception of the cross section differential in x 1 , for which the corrections change 
the cross section by up to 15%. 



3 Detector Description 

A detailed description of the HI detector can be found elsewhere [38] and only the components 
relevant for this analysis are described here. 

The HI central tracking system is mounted concentrically around the beam-line and covers 
polar angles 3 20° < 6 < 160°. The transverse momenta and charges of charged particles are 

3 The z axis of the right-handed coordinate system used by HI is defined to lie along the direction of the proton 
beam with the origin at the nominal ep interaction vertex. 
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measured by two coaxial cylindrical drift chambers [39]. Two drift chambers which provide 
accurate measurements of the z coordinate of charged tracks and two multi-wire proportional 
chambers which trigger on these tracks are placed on either side of the inner main drift chamber. 

The tracking system is surrounded by a finely segmented Liquid Argon Calorimeter [40], 
which covers the range of polar angles 4° < 9 < 154° and the full range in azimuth. It con- 
sists of an electromagnetic section with lead absorbers, 20-30 radiation lengths in depth, and a 
hadronic section with steel absorbers. The total depth of the calorimeter ranges from 4.5 to 8 
hadronic interaction lengths. The energy resolution obtained from test beam measurements [41] 
is a(E)/E « 0.11/ VE for electrons and a(E)/E « 0.5/ y/E for pions, with E in GeV. The 
absolute energy scale for hadrons is known for this analysis to a precision of 4%. A uniform 
axial magnetic field of 1.15 T is provided by a superconducting coil, which surrounds the calo- 
rimeter. 

The polar angle region 153° < 9 < 177.8° is covered by the SPACAL [42], a lead/scin- 
tillating fibre calorimeter with electromagnetic and hadronic sections. The energy resolution 
of the electromagnetic section is determined to be 0.07/v^E © 0.01 (E in GeV) [43]. Both 
calorimeter sections have a time resolution better than 1 ns. The SPACAL is used both to 
trigger on the scattered electron 4 and to measure its energy. In front of the SPACAL, an eight 
layer drift chamber, BDC [44], covers the polar angle region 151° < 9 < 177.5°. It is used to 
suppress background from neutral particles faking the scattered electron and, together with the 
vertex obtained from the central drift chambers, to measure the scattered electron polar angle 9. 



4 Data Samples and Event Selection 

The present analysis is based on data taken in the years 1999 and 2000, when electrons with an 
energy of 27.55 GeV collided with protons with an energy of 920 GeV. The data correspond to 
an integrated luminosity of 57 pb^ 1 . 83% of the data sample corresponds to e + p collisions, the 
remainder to e~p interactions. 

The kinematic region covered by the analysis is defined by cuts on the photon virtuality, Q 2 , 
the inelasticity, y, and by cuts on the hadronic final state. We require 

2 < Q 2 < 80 GeV 2 , 

0.1 < y < 0.85. 

The variables Q 2 and y were determined using the scattered electron energy and polar an- 
gle [15]. 

The final state has to contain at least two jets. Jets are found using the longitudinally in- 
variant k t jet algorithm [45] applied to hadronic final state "combined objects", boosted into 
the photon-proton centre-of-mass frame. The combined objects are constructed from tracks in 
the central track chambers and clusters in the SPACAL and LAr calorimeters in a procedure 
that avoids double counting [46]. The jet transverse energies, E? , and pseudorapidities, rj* , 

4 In the following, the notation "electrons" stands for both positrons and electrons. 
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are calculated relative to the 7*p collision axis in the 7*p centre-of-mass frame 5 . The jets are 
ordered according to their transverse energy, with jet 1 being the highest Ej, jet. 

The two jets with the highest transverse energies (leading jets) are required to have 

££i>7GeV, £* 2 >5GeV, 

-2.5 < r){ < 0, - 2.5 < rj* < 0. 
The asymmetric E^ cuts avoid regions of instability in the NLO calculations [47,48]. 

The reconstructed event vertex has to be within ±35 cm of the nominal interaction point, 
which substantially reduces contributions from beam induced background. To remove back- 
ground from photoproduction processes, a cut 45 < Ei(^j ~ Pz,i) < 75 GeV is applied, where 
the sum runs over all particles in the final state including the scattered electron. In total 105 658 
events satisfied the selection criteria. 

It is convenient to describe dijet events by means of the variable x^ ts , defined as 

7 " E (e*-p* z )' w 

hadrons 

where the sum in the numerator runs over the two leading jets and the sum in the denominator 
includes the full hadronic final state. Neglecting the masses of the partons and beam particles, 
the variable x^ ts represents a hadron level estimate of the fraction of the photon four-momentum 
carried by the parton involved in the hard scattering. 



5 Analysis Procedure 

The data were corrected for initial and final state QED radiation effects using samples of RAP- 
GAP events for direct photon interactions with and without QED radiation, processed through 
the full detector simulation and fulfilling all the cuts described in the previous section. The 
effects of trigger inefficiencies, limited detector acceptance and resolution were corrected for 
using an iterative Bayesian unfolding technique [49], which was applied to events generated by 
HERWIG and RAPGAP. For this purpose, 4 million events from each generator were passed 
through the full simulation of the HI detector and the same chain of reconstruction and analysis 
procedure as for the data. 

The binning for the final results was chosen such that the bin width is always larger than the 
resolution of the given quantity. The iterative Bayesian procedure converged in all bins of the 
measured quantities [21]. After unfolding, the correlations between neighbouring bins in the 
unfolded distributions were always less than 60%. The correction factors from the unfolding 

5 The pseudorapidity is defined by 77* = — ln(tan0*/2), where 9* is the polar angle of the jet axis with respect 
to the -y*p collision axis. Negative values of rf correspond to the photon fragmentation region. The pseudorapidity 
in the photon-proton centre-of-mass frame is shifted on average by -2.3 units with respect to the pseudorapidity in 
the laboratory frame. 
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procedure were cross-checked with a bin-by-bin correction method, performed using the same 
simulated MC samples, and agreement was found to within 5%. 

In all of the distributions studied, the presented cross sections are taken as averages of the 
cross sections obtained when correcting for detector effects using HERWIG and RAPGAP, 
since the description from the two models of the uncorrected distributions are of similar quality. 

The background in the event sample from photoproduction events, in which a hadron in 
the SPACAL is misidentified as the electron candidate, was estimated using PYTHIA [28] and 
PHOJET [50] MC samples of photoproduction events. This background is negligible for most 
of the bins and reaches 4% at the highest y. 

The systematic errors are added in quadrature. They are listed below in order of their size: 

- Model Dependence. The systematic error from the model dependence of the acceptance 
corrections is taken as half of the difference between the results when unfolded with 
RAPGAP and with HERWIG. This leads to an error of 5-10% on average, reaching 20% 
in the most extreme case. 

- Energy Calibration of the Calorimeters. Varying the overall hadronic energy scale of 
the LAr calorimeter by 4%, the hadronic energy scale of the SPACAL by 7% and the 
electromagnetic energy scale of the SPACAL by 1% leads to systematic shifts of the 
results by typically 10%, 2% and 4%, respectively. 

- Scattered Electron Angle. The polar angle of the scattered electron is measured with 
a precision of 1 mrad, which leads to a 3% (1%) systematic uncertainty in the lowest 
(highest) Q 2 region. 

- Trigger Efficiency. The uncertainty in the trigger efficiency leads to a 3% uncertainty in 
the measurement. 

- QED Radiative Corrections. 2% is taken as the systematic error in all bins [21]. 

- Stability of the Bayesian Unfolding Procedure. By varying the number of iterations used 
in the unfolding procedure, the uncertainty due to the unfolding instability is estimated to 
be typically less than 2% and at most 5%. 

- Photoproduction Background. The photoproduction background is subtracted statistically 
and half of the subtracted background is taken as the systematic uncertainty. 

- Precision of the Luminosity Measurement. The normalisation uncertainty due to the lu- 
minosity measurement is 1.5%. 



6 Results 

Differential cross sections for the kinematic region defined in Section 4 are discussed in the 
following sections and presented in Tables 2-5. 
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6.1 Comparison with NLO Calculations 



The triple differential dijet cross section is presented as a function of x^ ts in different bins of Q 2 
and E?p in Fig. 2. The variable denotes the transverse energies of the jets with the highest and 
second highest measured in the photon-proton centre-of-mass frame, so that each event con- 
tributes twice to the distributions, not necessarily in the same bin. The data are compared with 
the NLO direct photon calculations 6 performed with DISENT and JETVIP. The uncertainties 
from variations of the factorisation and renormalisation scales in the interval /x/2 to 2/x, as well 
as from hadronisation corrections, are illustrated. The scale uncertainties are typically around 
20%. Those from hadronisation are at the 7% level. We have also investigated the uncertainties 
due to variations of the proton PDF, using the prescription of [51]. The typical uncertainties are 
below 4% (not shown). Figure 2 demonstrates that the NLO direct photon calculations describe 
the data in the region of high x^ ts , where direct photon interactions dominate. For x J 7 cts < 0.75, 
the description is nowhere perfect, indicating the need for orders beyond NLO. The description 
of the data for x J 7 cts < 0.75 gets worse as Q 2 and decrease. The discrepancy is particularly 
pronounced at small x^ ts , low Q 2 and low E^ , where the data lie significantly above the theo- 
retical predictions, even taking into account the sizable scale uncertainty. The relative decrease 
of the cross section at low x J 7 ets as increases is of kinematic origin, due to the restrictions 
in the available phase space. Note that for x J 7 cts < 0.75, the JETVIP results are systematically 
lower than those of DISENT, whereas for x^ ts > 0.75 the opposite effect is observed. The 
discrepancy between DISENT and JETVIP is observed only for multi-differential distributions 
which include a jet variable. It gets substantially smaller for the inclusive dijet cross section 
d 2 a ep /dQ 2 dy [21] (not shown) and agrees within 2% for the total dijet cross section in our 
kinematic region. A similar level of agreement between JETVIP and DISENT was reported 
in [52] for the total dijet cross section. 

The data were also analysed in terms of jet pseudorapidities. Figure 3 presents the dijet cross 
section as a function of rj* in different bins of Q 2 and y, where rj* denotes the pseudorapidities 
of the jets with the highest and second highest E^ in the photon-proton centre-of-mass system, 
such that each event enters the distributions twice. The excess of the data over the theory at low 
Q 2 and low x J 7 cts observed in Fig. 2 is reflected in Fig. 3 in a similar excess at low Q 2 and high 
y, which is especially pronounced in the forward region of the laboratory frame (rj* ~ 0). 

Figure 4 shows the triple differential dijet cross section as a function of E^ in different bins 
of Q 2 and r/* , each event entering the distributions twice. The predictions of the NLO direct 
calculations agree well with the data at large Q 2 or at large Ej, for all r/* . On the other hand, the 
predictions clearly fail to describe the data in the forward region at low Q 2 and low Ej, . The low 
Ej, region is better described as rj* is reduced or Q 2 is increased. A similar discrepancy between 
the data and the NLO prediction has recently been reported for inclusive jet cross sections in a 
similar kinematic region [4]. The measurement in [4] indicated that the region where the NLO 
calculations fail to describe the data corresponds to the region where the ratio of NLO to LO 
predictions is largest. The same is true for the dijet cross sections. The NLO corrections are 
smallest in the backward region at the largest Q 2 and E^, where the ratio is approximately 1.1. 
The data are well described by the NLO direct calculations in this kinematic region, as can be 
seen in Fig. 4. On the other hand, the ratio of NLO to LO predictions for the forward region 

6 The resolved photon prediction of JETVIP, also shown in Fig. 2, is discussed in Section 6.2. 
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at small Q 2 and becomes as large as 9 and the data there are not reproduced by the NLO 
calculations. Corrections beyond NLO are therefore expected to improve the description in this 
region. 

The above comparisons show that in the region of low Q 2 , high y, forward r/* and low 
the data lie significantly above NLO QCD calculations for direct photons. This excess cannot be 
accommodated within standard theoretical uncertainties from scale variations and hadronisation 
corrections. 



6.2 Resolved Virtual Photons 

The pattern of the observed discrepancy between the data and the NLO calculations in Fig. 2-4 
suggests an explanation in terms of the interactions of resolved virtual photons, understood as 
an approximation to contributions beyond NLO. Of the NLO parton level calculations, only 
JETVIP includes a resolved virtual photon contribution 7 . Unfortunately, the dependence on the 
slicing parameter y c of the NLO JETVIP calculations of the resolved 7^ contribution (i.e. up to 
order aa z s ) is much larger than for the direct component [21], and the resulting calculations are 
therefore less reliable. In the absence of other calculations of this kind, the data are compared 
with the results of the full JETVIP calculations in Fig. 2 using y c = 0.003 (see Section 2.3), in 
order to see the qualitative effects of resolved photon interactions at NLO. 

The inclusion of a resolved 7^ contribution brings the NLO calculations closer to the data, 
though there is still a discrepancy between the data and calculations at low to moderate x J 7 ets 
and low Q 2 . The dominant part of the difference between the full NLO JETVIP results and 
the direct component comes from the O(aa^) term in the resolved 7^ contribution (Fig. 1 e). 
Including only the leading resolved 7^ terms (Fig. 1 c,d) has only a small effect [21]. 

6.3 Comparison with DGLAP Monte Carlo Models 

Unlike the parton level calculations, DISENT and JETVIP, discussed in the context of Fig. 2-4, 
LO MC models take initial and final state QCD parton showers into account. Their importance, 
together with the QCD improvements of (see Eq. (2)) and the simulation of soft under- 
lying interactions and hadronisation, is demonstrated in Fig. 5. The full HERWIG simulation, 
as described in Section 2.3, is compared with the HERWIG prediction without hadronisation 
or soft underlying event effects, and with the HERWIG calculation at the parton level without 
parton showers and with only the QED PDFs of virtual photons. The cross sections predicted 
by the full HERWIG simulation are in good agreement with the data in the low x^ ts region. The 
highest x^ ts region is not described so well. The largest difference between the cross section 
predicted by the full HERWIG simulation and that obtained from the parton level calculation 
comes from the initial and final state QCD parton showers, which effectively introduce an in- 
trinsic k T of the partons in the incoming proton. These effects increase the total dijet cross 
section in our kinematic region typically by 30% and by as much as 100% at low Q 2 , low 
and low x^ ts . Another 10% increase of the total dijet cross section arises from the change from 

7 Only the contribution of transversely polarised resolved photons is implemented in JETVIP. 
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QED to QCD-improved -Dj/ 7 *. Soft underlying events increase the total dijet cross section by 
4%. Their influence is largest in the region of low Q 2 , low E^ and low x^ ts , where the cross 
section is increased by 10%. A similar effect is observed when using the multiple interaction 
model implemented in PYTHIA. 

Previous analyses of jet production in low Q 2 ep collisions have compared with resolved 
virtual photon models that neglect longitudinally polarised photons. Figure 6 shows predictions 
of the direct and 7^, and 7^ resolved photon components separately. At high Q 2 the HERWIG 
direct photon prediction alone reasonably describes the shape of the a: J 7 cts distribution of the 
data, while at low Q 2 the resolved photon contributions are clearly needed. The contribution of 
longitudinally polarised resolved photons improves the agreement with the data. Not only do 
they increase the magnitude of the HERWIG predictions such that they become closer to data, 
but they also correctly reproduce the Q 2 and dependence. For a given interval of E^ , the 
ratio of 7^ to 7^ contributions increases with Q 2 , whereas keeping Q 2 fixed it decreases with 
increasing E^ . This behaviour is expected from Eqs. (5) and (6). Enhancing the PDF of 7^ in 
the resolved photon contribution by a constant factor does not lead to a comparably successful 
description of the data. 

As a result of the different y dependences of the photon fluxes in Eqs. (3) and (4), the dijet 
cross section as a function of y is different for longitudinal and transverse photons. Figure 7 
shows the event cross section as a function of y in different bins of Q 2 and x J 7 cts . In contrast 
to all previous dijet cross sections, each event contributes only once to the event cross section 
d 3 u ep /d(3 2 da; J 7 ets d?/ . The ratio of longitudinally to transversely polarised photons decreases 
with increasing y, as expected from the fluxes (see Eqs. (3) and (4)). The addition of the 
resolved longitudinal photon contribution brings the HERWIG predictions 8 closer to the data. 
The small contribution of 7^ compared to 7^ at large x^ ts is a consequence of the different x 7 
dependences of A/7* and A/ 7 £ (see Eqs. (5) and (6)). 

Figure 8 compares the measured dijet cross section as a function of r/* in different bins 
of Q 2 and y, presented already in Fig. 3, with the HERWIG prediction. The data are well 
reproduced by the complete LO MC model in shape. However, the absolute normalisation is 
not satisfactory, especially at low y. In agreement with the conclusion of Figs. 6 and 7, the 
importance of the resolved photon contributions increases in the forward jet region (rj* ~ 0), 
for low Q 2 and at high y. 

6.4 Comparison with CCFM Monte Carlo Model 

In Figs. 6-8, the data are also compared with the predictions of the CASCADE MC, employing 
the unintegrated PDF (set 1 of [36]) with the CCFM evolution equations. 

The CASCADE prediction describes the main qualitative trends in the data, except the Q 2 
dependence in the lowest E^ bin (Fig. 6) or at low x J 7 ets (Fig. 7). CASCADE also overestimates 
the data in the lowest y bin at high x^ ts (Fig. 7). On the other hand, CASCADE predicts a 

8 The low HERWIG prediction for all bins with x J 7 cts > 0.75 and at the lowest y for x J 7 0ts < 0.75 is partially due 
to a cut-off procedure in HERWIG, which suppresses the PDF of the virtual photon at large x J 7 ts . The resolved 7^ 
contribution of RAPGAP (not shown) leads to a rise with decreasing y that is similar to that in the data for the low 
£ J 7 ts range, though RAPGAP also lies below the data in the large x J 7 ts range [21]. 
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significant dijet cross section at low x^ ts (Fig. 6), much higher and closer to the data than the 
LO and NLO DGLAP predictions without the resolved photon interactions. Also, except for 
the highest Q 2 bin, dijet production in the forward region is reproduced better by CASCADE 
(Fig. 8) than by NLO direct photon calculations (Fig. 3). 

Large sensitivity of the CASCADE predictions to the choice of unintegrated proton PDF is 
observed [21]. The results shown here are based on set 1 of [36], where only the singular terms 
in the gluon splitting function are included. Switching to set 2 of [36], for which the full gluon 
splitting function is used, results in a reduction in the predicted cross section by up to 30 %. 
Set 2 gives the best description for different observables in another recent dijet measurement 
covering a similar kinematic region [3]. 



7 Summary 

Triple differential dijet cross sections in e ± p interactions are measured in the region of photon 
virtualities 2 < Q 2 < 80 GeV 2 and over a wide range of inelasticities 0.1 < y < 0.85. The 
data, covering the kinematic range E^, 1 > 7 GeV, ££ 2 > 5 GeV and pseudorapidities —2.5 < 
r)i, r)2 < 0, are compared with NLO and LO calculations, with and without resolved photon 
contributions or parton showers, as well as with a calculation based on k T factorisation and an 
unintegrated PDF of the proton. 

A sizable and systematic excess of the data over the NLO calculations of DISENT, which 
do not include a resolved virtual photon contribution, is observed for Q 2 < 10 GeV 2 , small 
jet transverse energies, E?, and small x^ ts , or equivalently, large jet pseudorapidities, rf. The 
excess observed for a; J 7 ets < 0.75 decreases with increasing Q 2 . 

NLO QCD calculations incorporating a resolved virtual photon, as implemented in JETVIP, 
bring the QCD predictions closer to the data, though there is still a deficit at low x^ ts , especially 
for low Q 2 . Unfortunately the JETVIP prediction for the resolved part of the dijet cross section 
is sensitive to the choice of the slicing parameter y c and must therefore be taken with caution. 

The significant role of initial and final state QCD parton showers, which are not taken into 
account in the NLO QCD calculations, is demonstrated. The inclusion of QCD parton showers 
in the HERWIG LO Monte Carlo model leads to a considerable improvement in the description, 
though a discrepancy remains in the region of high x J 7 ets . The best agreement with the data 
is obtained when both transversely and longitudinally polarised resolved virtual photons are 
included. 

CASCADE, which is based on the CCFM evolution scheme, does not involve the concept 
of virtual photon structure. The CASCADE description of the data is best in the region of 
moderate Q 2 between 10 and 25 GeV 2 . The Q 2 dependence of the cross section is less steep 
than in the data. 

To conclude, the data show clear evidence for effects that go beyond the fixed-order NLO 
QCD calculations. The importance of QCD parton showers and of the resolved contribution 
is demonstrated. 
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Figure 2: Triple differential dijet cross section, d 3 o~2j C t /dQ 2 dE?p dx J 7 ets , with asymmetric 
cuts (see text). The inner error bars on the data points show the statistical error, the outer error 
bars show the quadratic sum of systematic and statistical errors. Also shown are NLO direct 
photon calculations using DISENT (hatched area) and JETVIP (full line), as well as the sum of 
NLO direct and NLO resolved photon contributions of JETVIP (dashed line). All calculations 
are corrected for hadronisation effects. The inner hatched area illustrates the uncertainty due to 
the hadronisation corrections, the outer hatched area shows the quadratic sum of the errors from 
hadronisation and the scale uncertainty (shown only for DISENT). The scale factors applied to 
the cross sections are given. 
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Figure 3: Triple differential dijet cross section, d 3 a2j e t/dQ 2 dydr]* . Negative values of r]* cor- 
respond to the photon fragmentation region. See the caption of Fig. 2 for further details. 
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Figure 4: Triple differential dijet cross section, d 3 (T 2 jct/dQ 2 dr7* dE^. . Negative values of 77* 
correspond to the photon fragmentation region. See the caption of Fig. 2 for further details. 
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• HI data 
— HER dir+res T +res L (full) 
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Figure 5: The triple differential dijet cross section, d 3 a2jet/dQ 2 dE^ dx J 7 ets , from HI data, com- 
pared with the predictions of the full HERWIG simulation as defined in Section 2.3 (full line), 
HERWIG without hadronisation or soft underlying event effects (dotted line) and HERWIG at 
the parton level (dashed line). In the latter case, the QED PDFs of and (see Eqs. (5) 
and (7)) are used in the resolved photon contributions. 
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Figure 6: The triple differential dijet cross section, d 3 <72j e t/dQ 2 d.E^ dx J 7 ets , from HI data com- 
pared with the predictions of HERWIG and CASCADE. The dark-filled histograms show the 
direct HERWIG contribution. The light-filled histograms the resolved 7^ HERWIG prediction 
and the full line is the sum of all direct, 7^ and 7^ resolved HERWIG contributions. 
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Figure 7: Triple differential event cross section, d 3 a ep /dQ 2 dx^ ts dy . See the caption of Fig. 6 
for further details. 



26 



• HI data 



HERWIGdir — HERWIG dir+res T +res L 

... CASCADE 



0.1 < y < 0.25 



1 1 1 ^ i\ yri\j tin 

HERWIG dir+res T 

0.25 < y < 0.5 



0.5 < y < 0.85 




Figure 8: Triple differential dijet cross section, d 3 <T 2 j e t/dQ 2 d|/dr/*. See the caption of Fig. 6 for 
further details. 
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Table 2: Triple differential dijet cross section, d 3 <7 2 j e t / dQ 2 dE^ dx J 7 ets . The cross section is given 
together with the statistical and systematic uncertainties. The correction factors for hadronisa- 
tion effects applied to the NLO QCD predictions are also given. 
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Table 3: Triple differential dijet cross section, d 3 cr 2 j e t/dQ 2 d?/d?7* . See the caption of Table 2 
for further details. 
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0.301 


0.004 


0.041 


0.98 






7-10 


0.477 


0.005 


0.064 


1.03 






10-15 


0.238 


0.003 


0.030 


1.05 






15-20 


0.090 


0.002 


0.010 


1.02 






20-30 


0.0229 


0.0007 


0.0042 


1.00 



Table 4: Triple differential dijet cross section, d 3 o- 2 j C t/d<5 2 d?7* dE?. See the caption of Table 2 
for further details. 
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Q z 


x 7 


y 


d a cr C p/dQ 2 da;J 7 0ts dy 
(pb/GeV 2 ) 


^stat 


^syst 


hadr. corr. 


(GeV 2 ) 






(pb/GeV 2 ) 


(pb/GeV 2 ) 




2.0-4.4 


- 0.75 


0.10-0.25 


586 


12 


100 


0.98 






0.25 - 0.40 


409 


9 


63 


0.93 






0.40-0.55 


283 


7 


41 


0.88 






0.55-0.70 


199 


5 


27 


0.86 






0.70-0.85 


192 


8 


27 


0.87 




0.75 - 1 


0.10-0.25 


1360 


31 


200 


0.96 






0.25 - 0.40 


653 


15 


63 


1.11 






0.40-0.55 


289 


9 


25 


1.15 






0.55-0.70 


164 


7 


14 


1.13 






0.70-0.85 


89.6 


6.2 


7.3 


1.12 


4.4-10 


- 0.75 


0.10-0.25 


198 


3 


27 


0.99 






0.25 - 0.40 


126 


2 


17 


0.95 






0.40-0.55 


92 


2 


12 


0.90 






0.55-0.70 


66.9 


2.1 


9.1 


0.88 






0.70-0.85 


53.6 


2.2 


7.1 


0.88 




0.75 - 1 


0.10-0.25 


744 


11 


74 


0.97 






0.25 - 0.40 


291 


6 


30 


1.12 






0.40-0.55 


138 


4 


12 


1.17 






0.55-0.70 


69.0 


2.8 


5.5 


1.15 






0.70-0.85 


34.3 


2.4 


2.7 


1.15 


10-25 


- 0.75 


0.10-0.25 


58.8 


0.9 


8.6 


1.04 






0.25 - 0.40 


37.8 


0.7 


4.6 


0.97 






0.40-0.55 


25.0 


0.6 


3.1 


0.91 






0.55-0.70 


17.5 


0.5 


2.2 


0.91 






0.70-0.85 


13.4 


0.7 


1.5 


0.88 




0.75 - 1 


0.10-0.25 


294 


4 


30 


0.97 






0.25 - 0.40 


112 


2 


12 


1.13 






0.40-0.55 


49.6 


1.4 


4.0 


1.17 






0.55-0.70 


24.0 


1.0 


1.8 


1.14 






0.70-0.85 


12.35 


0.84 


0.91 


1.14 


25-80 


- 0.75 


0.10-0.25 


9.5 


0.2 


1.8 


1.07 






0.25 - 0.40 


6.98 


0.15 


0.85 


1.04 






0.40-0.55 


5.00 


0.15 


0.51 


0.97 






0.55-0.70 


3.32 


0.17 


0.36 


0.93 






0.70-0.85 


2.25 


0.27 


0.22 


0.92 




0.75 - 1 


0.10-0.25 


72.5 


0.9 


8.6 


0.98 






0.25 - 0.40 


31.4 


0.6 


2.9 


1.12 






0.40-0.55 


15.0 


0.5 


1.3 


1.14 






0.55-0.70 


8.99 


0.47 


0.88 


1.16 






0.70-0.85 


4.01 


0.57 


0.85 


1.13 



Table 5: Triple differential event cross section, d 3 a cp /dQ 2 dx^ ts dy . See the caption of Table 2 
for further details. 
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